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Abstract: We further analyze the definition and the calculation of the heavy quark 
impact factor at next-to-leading (NL) logs level, and we provide its analytical ex- 
^ . pression in a previously proposed fc-factorization scheme. Our results indicate that 

fc-factorization holds at NL level with a properly chosen energy scale, and with the 
^ . same gluonic Green's function previously found in the massless probe case. 
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1. Introduction 



Recent improvements [||] of the next-to- leading logo; (NLx) results in the BFKL 
framework, have stabilized the small-x behaviour in QCD, so that a phenomenological 
analysis of deep inelastic processes (DIS) seems now possible. 

However, both the gluon density (satisfying the improved equation) and the 
impact factors are needed in order to use fc-factorization (Sec. H) to compute DIS or 
double DIS processes. So far, NLx impact factors have been found for the unphysical 
case of massless initial quarks and gluons only Partial features for massive 

quarks P and for colourless sources P are known too. 

In this paper we derive complete results for the case of initial massive quarks, 
with a twofold purpose. First, we want to check the validity of the fc-factorization 
scheme introduced in Ref. [Q, or, in other words, to derive the same gluon Green's 
function with an explicit massive quark impact factor which satisfies the expected 
collinear properties. Secondly, we develop as a byproduct some analytical techniques 
which are needed to deal with two-scale problems, which are hopefully useful to cope 
with the physical cases also. 

The results of the paper rest on two observations. The first one, motivated 
in Sec. ^, is that the factorized scale relevant in a high-energy two-scale process 
coupled to heavy quarks is sq = Max{ki,mi)Max{k2,m2), where ki and /c2 denote 
the relevant gluon virtualities, rather than sq = kik2, as in the massless quark case. 



In fact, by subtracting the kernel contribution with such a scale we are able in Sec. ^ 



to derive a result for the massive quark impact factor which is finite for s ^ oo, and 
has all the desired properties. 

The second observation is that we are able to disentangle the (m//c)-dependence 
of the impact factor by explicitly computing its Mellin transform and its inverse. 
Given the singular energy dependence of the squared amplitude and of the phase 
space in the intermediate steps, this is by no means a trivial result and requires 
a careful handling of Mellin transform integrals in dimensional regularization, as 
explained in Sec. |0 . 

The outcome of such analysis is that the NLx constant iJ-kernels, previously 
introduced in the gluon Green's function [Q, are indeed probe- independent, and 
that the ensuing impact factors only contain factorizable single logarithmic collinear 
singularities. The use of such information in the improved small-x equation and the 
left-over problems are discussed in Sec. ||. 



2. k- Factorization in dijet production 

We consider the high-energy scattering of two partons a and b with momenta pi 
and p2 respectively. Following 0, the colour averaged differential cross section is 
factorized in a gauge-invariant way into a Green's function Q^^ and impact factors h^, 



1 



and K (Fig. |T]) 

do-ab 



d[fci] d[k2] 



duj 



(2.1) 



We adopt d[fc] = d?^'^^k/Ti^^'^ as transverse space measure. The transverse plane is 
defined with respect to the incoming momenta pi and p2- The transverse momenta ki 
and ^2 play the role of hard scales of the process. By definition, the impact factors 
are free of high-energy gluon exchanges, which are subtracted out, but can still 
contain collinear singularities which need to be factored out. The Green's function 
incorporates all the Regge-gluon exchanges between the two partons. The energy- 
scale So(fci, ^2) will be chosen later on. 

At the next-to-leading logx (NLx) accuracy the Green's function has the 
following general form 



'l + asHL) 



a. 



n -1 



1-—{K^ + Kml) 



'l + asHn) 



[2.2) 



where Kq and K^l are the leading logx 
(Lx) and the NLx BFKL kernels respec- 
tively, Hr{Hl) are operator factors introduced 
in HI and 



a. 



^?'r(l-£)/i 



2e 



TT 



(4vr; 



l+£ 



is the dimensionless strong coupling constant. 
As explained in H, the identification of 



the second order impact factors, h^^ and h^^ , 
is affected by a double factorization scheme 
ambiguity, due to both the choice of the scale 
So and of the kernels Hr{Hi). The latter were 
introduced by Ciafaloni and Colferai (CC) in |^ 
so as to provide partonic impact factors free of 
double log collinear divergences for the factor- 

ized scale choice Sq = kik2. It was also shown that the left-over single logarithmic 
divergences could be factorized by the usual DGLAP approach. A different factoriza- 
tion scheme, allowing double logarithmic divergences, was used instead in , where 
an integral representation for the massive quark impact factor was presented also. 

In this paper we extend the CC scheme to the massive quark case, by showing 
that the gluonic Green's function stays the same and that the collinear divergences 
of the massive impact factor stay single logarithmic too. Although the use of the H 



(1) 




Figure 1: Diagrammatic representa- 
tion of fc-factorization. 



2 



kernels is optional for colourless sources 0, for which they can be incorporated in 
the impact factors, we think that they help in stabilizing the coUinear behaviour of 
the gluonic Green's function, as already noticed in 0. 

In the following we use the notation pi,p2 (Ai, A2) for the initial parton's mo- 
menta (helicities) and the indices 3, 4 (possibly 5) for the final ones, with the Sudakov 
parametrization 

ki=Pi-P3 = zipi - p2 + ki± , 

(1 - zi)s 

k2= P2-P4 = -T, ^—rP^ + ^2P2 + k2± , 

(1 - Z2)S 

where we have introduced Sudakov variables Zi and transverse spacelike vectors ki± 
perpendicular to the plane of the initial particle momenta light-cone basis (pi,P2) 

2 2 

_ mf _ mi _ 

Pi = Pi H P2 , P2 = P2 H Pi , 

s s 

Pi-kj± = 0, p^ = 0, pf = m^, with D — 2 Euclidean components ki : kf = —kfj_ > 0. 
We also define q = ki + k2 as the transverse momentum of parton 5. For simplicity, 
we use in the sequel q = \q\, ki = \ki\ and we consider parton b as massless, 1712 = 0, 
mi = m. 



3. Factorization scheme and calculational procedure 



Let's consider first the high-energy scat- 
tering of two partons a and b where a = q is a 
heavy quark of mass m with real emission of 
an extra gluon g that we assume in the heavy 
quark fragmentation region (Fig. In terms 
of invariants, we work in the kinematical re- 
gion S2 = {q + Pif ^ zis > Si = (ps + qY ~ 
g^/zi, so that zi > q/y/s is the fragmentation 
phase space boundary. The Born differential 
cross section in this high energy region can 
be calculated in a straightforward way (e.g., 
by eikonal coupling to the incoming parton b) , 




Figure 2: Real ghion emission in the 
fragmentation region of quark a. 
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and results to be 



da, 



qgb 



Nr. 



dZi d[fci] d[fc2] 
■Pgq (^1 5 £ 



m 



[q"^ + 7ii?zf\ [{q — Zik2Y + rn^zW 
zl{l-zi)[{2q-zik2)-k2f 



^2 [q^ + ra'^zf\' [{q — ^1^2)^ + m?'zf\' , 

(1 - zi) [q-{q~zik2) +m'^zl] 



[kl + m'^zfj [q"^ + rn^zl] [{q — 2:1^2)^ + nrfizl 
^1(1 - z^fk2-(2q - fc2) [fc2-(2g - k^) - z^kl] 
kl [kl + m?zl] ^ [q2 + m?zl] [{q - ^1^2)^ + m?zl] 



where 



1 



7^^{z,,e) = — [l + {l-zrf + ezl] , 
is related to the quark to gluon splitting function, 



m - 1 



r(i 



(3.1) 



(3.2) 



(3.3) 



k' /i2. 

is the leading order impact factor, being the same for quarks and gluons, is the 
renormalization scale, and 



A, = k^ h'^^^k) 



a. 



(3.4) 



r(i - 5)/i2- 

is a constant that contains the dependence on the strong coupling constant and some 
colour factors. 

Though complicated at first sight, eq. (PTT|) has some simple features that we now 
comment upon. First of all, only the Nc part is really relevant to our purposes, the 
Cp part being canceled with virtual correction upon zi and ki integration (see the 
following). We consider 

dcTqgb 



d^i d[A;i] d[fc2] 

■Pgq (^1 ; ^ 



= Ah^^\k2) 
(1 - zi) [q-{q- ziki) +m^zl] 



m 



[kf + w?zl\ [q^ + rn^zl] [{q — ^1^2)^ + rn^zl 
2i(l-Zi)2fc2-(2q-fc2) [fc2-(2qr-fc2) ^ 



ki [kl + m?zl] [q2 + m^zl] [{q - 21^2)^ + m^zl] 

The latter expression reduces, as expected to the known result for m 
matches the Lx differential cross section 



qgt> 



dZid[fci]d[fc2] 



hf{k{)ht\k 



a. 



1 



q^Vil - £)/i2^ zi 



(3.5) 
0, and 

(3.6) 
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in the limit Zi 0. However, for eg. ( P-6|) to be a good approximation to eg. ( |3.5|) , 
we should reguire 

Zi<^q/ki, Zi<^q/m, Zi<^ki/m. (3.7) 

The first two cutoffs can be summarized by Zi < q/Max{ki, m), which is a coherence 
condition for the case of heavy guarks, saying that the rapidity of the gluon cannot 
exceed that of the final guark. By integrating the leading expression (|3.6|) with the 
constraints ( |3.7| ) in the fragmentation region zi > qj^fs-, we obtain 

1 da^g _ ashf{k,) ( y/i q_ \ 

ht\k^) d[fci] d[fc2] q2r(l - e)y?- V Max(A;i,m) k^^'^'"') ' ^^'''^ 

This expression is an estimate of the leading contribution contained in the complete 
result eg.( p.l|) , which should be subtracted out in order to yield the impact factor in 
the massive guark case. 

Compared to the subtraction (or factorization) scheme adopted in for m = 0, 
the expression ( |3.8| ) differs by the replacement k^ — * Max{ki,m), which leads, by 
adding the symmetrical fragmentation region, to the choice for the factorized scale 
in eg. (pTTP 

So = Max{ki,mi)Max{k2,'m2) , (3.9) 

mi being the mass of guark a and m2 the mass of guark b. 

By considering now both real and virtual contributions to the fragmentation 
function Fq{zi, ki, ^2), we are led to introduce the following definition of the impact 
factor hq\k): 



/ dzi / d[ki]Fq{zi,ki,k2) 

(3.10) 

where 



a,Ko{k,M)= 2rn"^ ^ + 2^W(fc^)^[q] , 6[q] = -n'^^ 5'^''' {q) , 

(3.11) 

is the leading kernel, with 

(i).,2^ g''N,k^ f d[p] asT\l + e) fe^" 



(47r)2+- y p2(fc_p)2 - 2£r(l + 2£)U'/ ' ^^'^^^ 
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the gluon Regge trajectory. Eg. (|3.10| ) reduces for m = 0, to the definition adopted in 
Ref. 0, and in particular contains the subtraction term logq/kiQgk^ which provides 
the expression {Hr = H\ = H) 



(3.13) 



for the H kernel in the fe-factorization formula. 

In order to simplify our subsequent calculations, we shall then use the known 
result 01 for m = 



11 



^ 3 



1 

-e 
2 



67 
18 



5n/ 
9iV, 



(3.14) 



and we shall explicitly compute only the difference for a non vanishing mass. For 
this reason, we write the fragmentation vertex for a heavy quark as the massless 
fragmentation vertex plus an extra quark mass dependent contribution that cancels 
out for m = 



Fq(zi, fci, fca) = i^<r^°(^i, fci, ^2) + AFq(Zi, fci, 



(3.15) 



Then, we find the following relationship between the massless quark impact factor 
and the heavy quark impact factor 



^S'Hfc2) = ^S;Lo(fc2) + d^i j d[fci]AFq(^i, fci, k, 

+ j d[ki]ash^;\k,)Ko{k,,k2) log^e^fc, . 



(3.16) 



The most complicated integral that remains in the r.h.s. of eq.( p.l6|) is the one for 
AFq which will be calculated through its Mellin transform in k2. Notice also that 
the integration limits in zi have been extended down to zi = 0. Since AF^ is regular 
at 2;i = this change introduces only a negligible error of order 1/s. 



4. Mellin transform and its inverse 

In order to perform the calculation outlined in eq.( p.l6| ), we proceed in two steps. 
First, we perform analytically the ki integration of eq.( p.l| ) by reducing the ki- 
integrals to two denominators, as explained in detail in Appendix |^. Then, we 
consider the virtual contributions [0 quoted in Appendix 0, and we organize them 
in terms of momentum fraction integrals only. Finally, summing up real and virtual 
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contributions to the fragmentation vertex (see eq.( |A.6| ) and eq. (p.6| ) at the appen- 
dices) (hehcity non conserving not included) we obtain the following expression for 
the difference AFq(fe2), arising from the second term in the r.h.s. of eg. ( P .161) 



AFq(fc2) 



AK 



[m 



2\£ 



2(1 + 2e) kl 



+ 



m-e) 



1 rl 



^0 



dzi dx 



l + e 



Zl 



2 
1 



-Zl 



\x{l — x)k\ + im?z"'^^ ^ 



1 pI 



zi(l — Zl) dzi dx 



[x(l — x)kl + m?zl\ 



l~e 



(4.1) 



where in Arq^'*(fc2), eq.( |B.5D of Appendix 0, the integration variable x has been 
identified with zi to simplify the sum. Note again that, because of the subtraction 
of the m = part, the Zi-integrals are convergent at Zi = 0. 

4.1 Mellin integrals 



To perform the last integrations in eq. ([4.1|) we first calculate its Mellin transform 



AFq(7) = r(l + £) {m 
yielding 

AFq(7) = A, (m 



2\~e 



d[k2 



1,2 \ 7-1 



AFq(/e2) 



T{^ + e)T{l-^-2e)T\l-^-e) 



l + e 



+ 



8r(2 - 27 - 2e) 
1 1 



7 + 2£ 1 - 27 - 4e V 1 - 7 - 2e 3 - 27 - 2e 



(4.2) 



Although this Mellin transform is finite for £ — > the limit e = cannot be taken 
in this expression. In fact, it is straightforward, though not trivial, to show that the 
r.h.s. of eq.( |4.1| ) behaves as 



(^2)^ ^ , 



fc2<»Tt 



{kl)-\m^-^ 



fc2>m 



(4.3) 



Therefore, the Mellin transform converges only in the small band l—2e< Re'~f < 1—s 
and the e dependence should be kept until the end. 

To recover AFq(A;2) we should calculate the following inverse Mellin transform 



AFq(fc2) = 



l-2e<Re7<l-e 



d7 

27ri V rn? 



T^Azr^] AFq(7) . 
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We consider first the limit > m?. Then, we displace the integration contour 
around the positive real semiaxis, enclosing all the poles placed in 7 > 1 — e, the 
smaller one giving the smaller power in m/k2. The first pole, at 7 = 1 — e, yields the 
following result 

AF,(fc,) = a. hf[k,) I - (^)V 0{m/k,)^ . (4.4) 

For ^2 < we displace the integration contour around the negative real semiaxis, 
enclosing all the poles placed in 7 < 1 — 26. As before, the first pole, at 7 = 1 — 2e, 
gives the answer at order 0{k2/m) 

AF,(*.)^ftr(..,{.».(.^)[-i±^-.og(|) 

+ V(l -e)- ^(1) - 2^{e) + 2^{2e) + O^k^/m) | . (4.5) 

Since ui^^^ ~ 1/e is infrared singular (cf. eg. ( p.l2| )), both formulas show double 
logarithmic singularities of type and 1 / e\og{k\/ m?) . 

The last ingredient we need in order to extract the next-to-leading quark impact 
factor is the last term in the r.h.s. of eq.( |3.16| ). For the real emission contribution 
to Kq we get the integral 



'q^T{l-e)n^^ ^ki 
We use the following representation 

valid for a, 6 > 0, which allows us to write 



2 J J {k\Y+^ 

A, r(l + A - e)V{e)V{e - A) . 2,a,^2,-i-a+. 



2 J r(i + A)r(2;- A) . (4.7) 

The integrand vanishes for |A| — 00 in all directions apart from the real axis. As 
before, we consider first the case k\ > m? and displace the integration contour around 
the positive real semiaxis enclosing all the poles placed in A > 0. The smaller pole, 
at A = e, gives us the result at order 0{m/k2) 

I. - ^. ("4 2.^r(i - ')r(i + .) )' + ■ ^'■''> 
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On the other hand, for fcg < m^, we consider the poles placed at the negative real 
semiaxis and therefore the residue at A = —a with a — >■ 0"*". By taking into account 
also the virtual contribution to in this case, we obtain 



(4.9) 



+ 2[^(l)-^(l-£) + ^(£)-^(2£)] 

Eq.( [4.8|) and eg. ( [4 .91) show double log singularities also. 
4.2 Impact factors 

Finally, summing up all the pieces according to eq. ( p.l6|) , the impact factor for heavy 
quarks at the next-to-leading level can be written as 

where 



(4.10) 



,,K.a)(^2. r(i + 2.) 



2 + 3£ + 2^2 



2(1 + 2e) r(l-£)r(l + £) 



0{m/k2 



(4.11) 



is valid in the limit fcg > m?, and 



^(l)-^(l-£)- 



1 + 5e - 2e 



21 



+ 0{k2/m) 



(4.12) 



2(1 + 25) 

is valid for k\ < m^. 

We notice, in the first place, that all double log contributions of type 
and 1/e log(fc2/'ni^) appearing in eqs. (|4.4| - |4.5|) and (|4.8| - |4.9| ) have canceled out in 
eqs. ( |4.11|) and (|4.12| ). This means that indeed our subtraction of the leading kernel 
was effective, thus lending credit to the scale ( |3.9| ) and to the if-kernel ( |3.13| ). 

The remaining singularities of the impact factor are single logarithmic ones ~ 1/e, 
having the structure 



^S'nfc2)|_^ = h^^\k2)(^lu;^'\kl) - l.;«(m^)e,,^ - ^^«(/c^)e^.. 



(4.13) 



Here the first term has the customary collinear interpretation 0], as coming from 
the finite part of the q — g anomalous dimension 



7gq 



TTUJ 



CpCXs 

27r 



3 1 
2-r 



(4.14) 
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while the remaining ones - depending on the scale Min{k2,m) - do not have such 
interpretation. Note, however, that a finite mass scale change m ^ cm, will produce 
exactly this type of contributions from the singular integration of the Ko kernel 
in eq. p.lOp acting on hq\ki) over the region < ki < Min{k2,m), leading to the 
expression 



UJ 



m k2 



(4.15) 



Therefore, the singularities ( |4.13| ) can finally be interpreted in the form 



^S'Hfc2)|_^ = hf{k,) ^ [uj^'\kl) - ^^(m^)] e.,^ + 5h^\k,) 



I smg 



hf{k,) 



2 

asNc 
27r 



^log^)0fc,^ + 5/lS')(fc2) 



(4.16) 



meaning that h^^ is actually finite, with the log(A;2/m^) dependence predicted by 
the DGLAP equations, apart from a proper mass scale change in eq. (|3.13|) . In other 
words, the scale leading to a finite massive quark impact factor differs from eq. (|3.9| ) 
by a finite renormalization of the quark mass, which is a normal ambiguity in this 
type of problems. 

Our final result for the heavy quark impact factor at the next-to-leading level 
reads 



h,ik2) = h(^\k. 



I smg 



h„(kn 



"qV^iJl finite ' 

where the singular piece is defined in eq.( [4.16 ) and 

Mk2)\f,nUe = h^:\^sik2)) 



(4.17) 



xa + 



2-K 



^ Re7>l ^ 



2+ Res[h{i)]\Qmk2 I 



(4.18) 



Re-y< 

is the finite contribution to the heavy quark impact factor, with 
67 7r^ 5n/ 

18 ~ y 



(4.19) 



the constant term of the impact factor for massless quarks, eq. (|3.14| ), and 



Hi) = I 4 



m 



1-7 



(1-7)^ 



r(7)r^(i-7) 
4r(2 - 27) 

l-7)+V^(7)-2^(l)] 



1 2 

7 ^ 1 - 27 



7 3-27 



(4.20) 
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The compact expression ( [4.20I ) was obtained by adding to the Melhn transform (|47 
the contribution of eg. ([4 .71), with A — > — 1 + 7 + 2e, and the Melhn transform of 



the virtual piece of the last term at the r.h.s of eq. (|3.16|) . The sum of these three 
terms, apart from 7 = 1 — e, 1 — 25, whose contributions have already been treated 
separately, is finite and therefore was expanded for e ^ 0. As for the massless case, 
the singularities proportional to (11/6 — Uf/SNc), the beta function, were absorbed 
by the running strong coupling constant 0:5(^2). The function h{'~f) provides the 
corrections of order 0{m/k2} and 0{k2/m) to the impact factor for m? < k\ and 
> kl respectively, yielding the following final result 



Rey>l 
X 



klJ r2(n + l) 

^ ' n=l ^ ' 

11 1 6 



2{n + iy 2(2n-l)2 (2n + 1)3 2(2n + 1)^ 
2 11 3 3 



n n + 1 2(2n-l) (2n + 1)2 2(2n + 1) 
X (^i;{n + 1) - ^{2n + 2) - 1 log | , (4.21) 



and 



i?e7<l ^ ^ 




^ T{2n + 2) \ J \n n+1 2n + 1 2n + 3 J ' ^ ' ' 
5. Conclusions 

Starting from the explicit squared matrix element for gluon emission in eq.( p.l| ) we 
have motivated the subtraction of the leading term in eq. ( |3.8| ), and we have performed 
the kl and Zi integrals needed to provide an explicit result for the heavy quark impact 
factor in eq. ( [4.16| ) and eq. ( [4.18| ). 



Even if the cross section being investigated is unphysical, the relevance of our 
results stems from the consistency of the following features: (i) the validity of the 



fc-factorization formula ( p.l|) with scale Sq = Max{ki,mi)Max{k2,ni2)', (ii) the ex- 
plicit expression of the impact factor with factorizable single logarithmic coUinear 
divergences, and (iii) the probe-independence of the subleading if-kernels of the CC 
scheme defined in eq.( |3.13| ). 

Even if such universal extra kernels can be reabsorbed in the impact factors for 
colourless sources 0, they help clarifying the structure of the coUinear limits for 
two-scale processes as elaborated at length by CCS pi]. Here, it was shown that the 
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gluonic Green's function including such kernels, is free of double logs of collinear type 
for both ki ^ k2 and A;2 3> fci. As a consequence, even colourless impact factors 
will show, in the present scheme, simple collinear properties, as expected from their 
DGLAP analysis [g. 

Of course, the real problem is to provide an explicit expression for the DIS 
impact factors. But - if the lesson learned form the L and NL kernels is still valid 
- the impact factor's magnitude is not expected to be much different from their 
approximate collinear evaluation. 



Aknowledgements 

We wish to thank Gianni Camici for very valuable contributions in the early stages 
of this paper. Work supported in part by the E.U. QCDNET contract FMRX-CT98- 
0194 and by MURST (Italy). G.R. acknowledges financial support from INFN (Italy) 
and BMBF Project 05HT9VKB0 (Germany). 



A. Real contribution to the fragmentation vertex 



The differential cross section for real gluon emission off a heavy quark, eq.( p.l|) , can 
be simplified by using identities of the following type 



^3(1 



Zi) [{2q - Zik2)-k2\ 



[q^ + m?zf\' [{q — 21^2)^ + rTi?z\f' 
1 

+ 



2:1(1 - Zi] 
2 



1 



2.212 



[q^ + m'^zi 



[{q - z^k^Y + w?zlf [q" + m?zl] [{q - Zik2Y + w?zl] 



(A.l) 



to split the full expression into several contributions with at most two different prop- 
agators free of q dependences at the numerator. After some algebra, we obtain the 
following expression for the real contribution to the fragmentation vertex of quark q 



+ 



A. 



m 

^i(l-^i)-7j 



1 



2.212 



[{q - Zik2y + m?zl 
2 



[k\ + im?zl\' 



2/ [q^ + 7ii?zf\ [{q — Zik2Y + im?zf\ 



+ 



ygq(zi,g) + 2zi(l-Zi)(mVfc^) 

2[kl + m^zl] [q^ + m^zj] 
(1 - 2i)^ygq(^i,e) + 2^1(1 - ^i)(mVfc^) 

2 [kl + m^zf] [{q - ^1^2)2 + m^zf] 



(A.2) 
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According to eq.( |3.10| ), this expression has to be integrated first in ki and then in 
zi, for zi > q/y/s. To perform the first integration we use 



/ 



d[y 

[kl + m^zfY [{ki + pf + m?zl]'^ 



V{a + [3-l-e) 

r(a)r(/3) 



dx 



(1 — X 



1/9-1 



(A.3) 



then, we obtain 



F^,real{zi,k2) = j <\[ki]F^^real{zi, ki, k 

= Aj:{l-e 



2\e 



ki-) 



+ l\x^ zj'?^,izue) + 2z^'+'%l-z,)imyk 







2 



[x(l — x)kl + m^] 



l-e 



+ 



+ 



ygq(^i,e) + 2^i(l-^i)(mVfc^) 

2 [x{l-x)kl + m^zlY~' 
(1 - z,y?^^izu£) + 2^i(l - ;2i)(mVfc^) 



2 - x)(l - + m'^zf] 



l-e 



(A.4) 



By subtracting the massless contribution 



AF^^realiZi, ^2) = F^^reali^i, ^2) - i^q^^azl^^l, ^2 



(A.5) 



we get an expression that is regular at = and therefore can be integrated down 
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to Zi = without changing the final result. Let's define 



Jo 



Cf l\ T(-e) (m 



2\e 



2 1 l + 2e kl 



1\ Ti-e) 



1 



l + 2e 2 



X 



dx 



[x{l — x)kl + m?'\ ^ ^ 
2(mVfc^) dx 



[x{l 



^)kf' 



l-e 



+ 



m-e) 



1 rl 



JO 



dzi dx 



zi l + e 
H —Zl 



2 
1 



+ 



\x{l — x)k\ + m^z^] 



[x(l-a;)fc2] 



21 1-^ 



[x(l-x)(l-Zi)2fc^+m2z2]l " -X)(l 



,21 l-e 



+ 



2 /•! /•! 



2 "'0 ^0 



+ 



2i(l — d2;i dx 
1 



— x)(l — ZiYk^ + m?zl\ 



l-e 



(A.6) 



Notice that some of the integrations has been kept undone. This long expression, 
although cumbersome, will be drastically simplified after adding the virtual contri- 
bution before doing any further integration. 



B. Virtual contribution to the fragmentation vertex 



The correction to the cross section due to virtual emission, including subleading 
effects, for general parton-parton scattering, can be extracted from the amplitude of 
Ref. i (t = -fc?) 



(-) 

aa 



1 1 
t 



1 + uj{-t) log 



-t 



r 



(+) 



5a 



bb PA4,A 



bb "A4,-A2 



(B.l) 



were r(+) and r(-) are the helicity conserving and the helicity non- conserving con- 
tributions respectively. At the order we are working, the virtual terms contribute to 
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the fragmentation vertex as follows 



(B.2) 



As for real emission, we split rq^"* into a massless contribution plus a quark mass 
dependent extra term 



r(+) 



qq 



qq,m=0 ' qq ' 



where AFq^'' cancels for m = 0, being 



.(+) 

qq,m=0 



(k) 



+ 



1 + 2e V4(3 + 2e) 
1 



2 



7 1 + e 



4 A'c 3 + 2£ 



1 
2 



2 Cf 
^ el^ \ l + 26 



e 
+ 2 



and 



ArW(fc) 



4r(l - 



X 



+ 



da; 
2(mVfc^ 



1 



l + 2e 



+ 



1 



— x)k'^ + m^] 



l-e 



dx 



Cf 2Vi-E) im^y 



1 + 2£ Jo [a;(l - + m2] 



l-e 



rri 



1 /•! 



JQ 



l + 2£ k^ 



X I + e 
1 — X 



X 



X 



+ 



l[y{l-y){l-xyk^ + m^x^Y ' [y{l - y){l - x^k^f ' 
2m^ x{l — x)dxdy 



'0 — y)(l — x)2fc^ + m^x^] 

We define 

AFq,„irt(fc2) = y dzi j d[ki]^F^.,irt{zi,ki,k2) 
= ;i(°)(fc2)2ArW(fc2) . 



l-e 



(B.3) 



(B.4) 



(B.5) 



(B.6) 



Adding up together the real, eq.( |A.6| ), and the virtual, eq. (|B.6|) , contributions to 
AFq(fc2) most of the terms cancel. To see this cancellation it is enough to identify 
the integration variable x appearing in Arq^''(fc2), eq.( [B.5| ), with the momentum 
fraction zi appearing in eq. ( |A.6| ). In particular, the Cf contribution fully cancels 
and the part gives the the simplified result presented in eq.( [4.1| ). 
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